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ABSTRACT 

A numerical study was performed to investi- 
gate the unsteady, multidimensional flow inside 
the combustion chambers of an idealized, two- 
dimensional, rotary engine under motored condi- 
tions . The numerical study was based on the 
time-dependent, two-dimensional, density- 
weighted, ensemble-averaged conservation equa- 
tions of mass, specie, momentum, and total 
energy valid for two-component ideal gas mix- 
tures. The ensemble -averaged conservation 
equations were closed by a K-c model of 
turbulence. This K-c model of turbulence 
was modified to account for some of the effects 
of compressibility, streamline curvature, 
low-Reynolds number, and preferential stress 
dissipation. Numerical solutions to the con- 
servation equations were obtained by the highly 
efficient implicit-f actored method of Beam and 
Warming. The grid system needed to obtain 
solutions were generated by an algebraic grid 
generation technique based on trans finite inter- 
polation. Results of the numerical study are 
presented in graphical form illustrating the 
flow patterns during intake, compression, 
gaseous fuel injection, expansion, and exhaust. 

DIFFICULTIES WITH THE SUPPLY of aviation gaso- 
line in several parts of the world have given 
rise to the need for more fuel-efficient and 
multifuel-capability general aviation aircraft 
engines (1-3).* One engine that is well suited 
for general aviation aircraft and can meet the 
requirements of higher fuel-efficiency and a 
greater tolerance to fuels with broader specifi- 
cations is the stratif ied-charge rotary engine 
(3). 


^Numbers in parentheses designate ref- 
erences at end of paper. 


In order to design more fuel-efficient and 
more fuel- tolerant stratif ied-charge rotary 
engines, it is necessary to have a good under- 
standing of the physics (fluid flow, heat tram; 
fer, and combustion) that take place inside the 
engines* combustion chambers. Relatively few 
investigators have studied the physics taking 
place inside rotary engines. Experimental 
research in rotary engines so far has been con- 
fined to studies on engine performance and 
emissions (4-10) and to high-speed photographic 
studies of combustion (11,12). Mathematical 
research in rotary engines has been broader in 
scope but has been limited to quasi- dimensional 
(13-22) and one- dimensional (23,25) analyses. 

To date, no studies have reported details of the 
unsteady, multidimensional flowfield inside 
rotary engine combustion chambers. 

The unsteady, multidimensional flowfield 
inside a rotary engine combustion chambers is 
complex. Owing to this complexity, progress 
towards understanding the details of this flow- 
field can only be made in a piecemeal manner by 
focusing on specific problems. Furthermore, it 
is necessary to first study simplified rotary 
engines and later study increasingly more 
realistic rotary engines. 

The objective of this study was to use a 
numerical method to investigate the flow pat- 
terns inside one of the combustion chambers of 
a simplified rotary engine under motored condi- 
tions. The simplified rotary engine studied 
here was two-dimensional (i.e., it was a rotary 
engine of infinite width) with "planar" gaseous 
fuel injection into the combustion chamber dur- 
ing the compression process. Also it was 
assumed that there was no leakage across the 
seals. The present study constitutes the first 
attempt at numerically studying the unsteady, 
multidimensional flow inside rotary engines. 

This study was conducted for a motored configu- 
ration in order to compare model results to 
experimental data which will soon be available. 
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In this report the simplified rotary engine 
studied is described in detail. Afterwards, 
formulation of the problem along with the 
essence of the numerical method of solution are 
presented. Finally, results for one numerical 
simulation are presented illustrating the flow 
patterns during intake, compression, gaseous 
fuel injection, expansion, and exhaust. 

DESCRIPTION OF PROBLEM 

A motored, two-dimensional rotary engine 
is shown in Fig. 1. It consists of a housing, 
an intake port, an exhaust port, a gaseous fuel 
injector, a rotor, three seals, three combustion 
chambers (the three regions bounded by the 
rotor, housing, and seals), and a shaft. The 
equations that describe the geometry of this 
rotary engine are given in Refs. 26 and 27 and 
summarized in Table X. The Cartesian coordi- 
nates X-l and Y 1 describing the inner surface 
of the housing (denoted as surface 1 in Fig. 1) 
are given by Eqs . (1) and (2) (Table I). The 
Cartesian coordinates X 2 and Y 2 describ- 
ing the outer surface of the rotor (denoted as 
surface 2, 3, and 4 in Fig. 1) are given by 
Eqs. (3) to (12) (Table I). The geometry of the 
rotor pockets is given by Eqs. (10) to (12). 
Here, it is noted that the parameter C rep- 
resenting the distance between the rotor apex 
and the housing is set at 0.004 m which is some- 
what large. The value chosen is to facilitate 
the numerical method of solution. Since leakage 
was not considered, the value chosen for C 
should not seriously affect the flow patterns. 
The origin and orientation of the Cartesian 
coordinate system is shown in Fig. 2(a). 

For this motored rotary engine, as the 
rotor rotates inside the housing, each of the 
three combustion chambers goes through five 
processes in the following order: intake, com- 

pression, gaseous fuel injection, expansion, and 
exhaust. Because of this symmetry, we only need 
to study the physics taking place inside one of 
the three combustion chambers. The combustion 
chamber chosen is the region bounded by surfaces 
1, 2, 5, and 6 in Fig. 1. 

As mentioned previously, leakage across the 
seals was not considered. If leakage across the 
seals was considered, we could still study the 
physics of the airflow taking place inside the 
rotary engine by only analyzing one combustion 
chamber . 

EXHAUST INTAKE 



Figure 1. - Geometry of a two-dimensional 
rotary engine. 



X-Y-t coordinate system. coordinate system 

Figure 2. - Spatial domain the the coordinate system. 


Initially, the combustion chamber to be 
studied is located at the beginning of the 
intake process (i.e., set 0 in Eqs. (3) and 
(4) to zero) and is filled with stagnant air at 
a constant static temperature T^ and pressure 
P£. Suddenly, the shaft starts to rotate at 
a constant angular velocity Q of 5000 rpm 
(revolutions/min) which in turn causes the rotor 
to rotate in the clockwise direction. This 
starts the intake process. 

During the intake process, air enters the 
combustion chamber through the intake port. 

The stagnation temperature T 0 and pressure 
Po (P 0 > p i» i.e., turbocharged) of the 
entering air are both constants. During the 
early part of the intake process, the exhaust 
port is also open. Thus some of the air also 
leaves the combustion chamber through the 
exhaust port. The back pressure of the exhaust 
port is maintained at P e xh ( p exh < p i) • 

Once the trailing apex (denoted by A^ in 
Fig. 1) passes point B on the housing (see 
Fig. 1), the intake port is closed. The reason 
for closing this port earlier is to facilitate 
the numerical method of solution. Once the 
intake port closes, the compression process 
starts. Later, between crank angles 0^ 
and 02 , a gaseous fuel (octane) is injected 
into the combustion chamber. Still later when 
the volume of the combustion chamber starts to 
increase, the expansion process starts. Once 
the leading apex (denoted by A L in Fig. 1) 
passes point C on the housing (see Fig. 1), the 
exhaust port opens initiating the exhaust 
process. The late opening of the exhaust port 
is to facilitate the numerical method of solu- 
tion. As noted earlier, the back pressure of 
the exhaust port is maintained constant. 

For the problem just described, we are 
interested in knowing the details of the flow 
patterns inside the combustion chamber during 
intake, compression, gaseous fuel injection, 
expansion, and exhaust. 

FORMULATION OF THE PROBLEM 

The equations that govern the problem des- 
cribed in the previous section are the conserva- 
tion equations of mass, species, momentum, and 
total energy for two-component ideal gas mix- 
tures. In this numerical study, the following 
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TABLE I. - EQUATIONS DESCRIBING GEOMETRY OF ROTARY ENGINE 3 


Equation Eq. No. 

X x » E cos (3A) + (R + C) cos (A) (1) 

Y 1 m E sin (3A) + (R + C) sin (A) (2) 

*2 - E sin C©) + X R0 cos (0/3) + Y ro sin (0/3) (3) 

Y 2 = E cos (0) + Y ro cos (0/3) - X R0 sin (0/3) (4) 

whore 


0 = 



Q dt = Qt 


(5) 


X R q = X R cos (t/6) + Y r sin (v/6) (6) 

y R 0 88 y R cos + Xr sin <*/*> (7) 

X R = R cos (2V) - (3E 2 /R) sin (6V) sin (2V) 

[ 1 1/2 

1 - (9E Z /R Z ) sin (3V) cos (3V) cos (2V) - P* (8) 

Yr = R sin (2V) - (3E 2 /R) sin (6V) cos (2V) 


+ 


2E 


[- 


2 2 2 
(9E /R ) sin 



cos (3V) cos (2V) - P 

y 


(9) 


P x = P cos 
Py - P sin 

/ 


p = < 


V 


(2V) 

(2V) 

0 



0 


*0^ - V)/^ 
XV - V 3 )/(V 2 

n(V - V 4 )/(V 3 



V, < V < V 
1 “ ~ T 

v 2 < v < v i 

V < V < V 
3 “ ~ 2 

V 4 * V * V 3 


( 10 ) 

(ID 


( 12 ) 


a A and V are parametric parameters. A varies between 0 and 2*. 

V varies between V L = *76 and V x = */2 for surface 2 (see Fig. 1); 
between = 5*/6 and V x = 7*/6 for surfaces 3; and between V L - 3*/2 
and V x = ll*/6 for surface 4. E = 0.015 m, R = 0.1045 m, C = 0.004 m, P^ = 

P 2 = 0.005 m, V X = V L + 0.75 (V T - V L ) , V 2 = V L + 0.625 (V x - V L ) , V 3 = V L = 

0.375 (V T - V L ), V 4 = V L + 0.25 (V x - V L ) , and Q = 5000 rpm. 


simplifications were made to the conservation 
equations: (I) all thermodynamic and transport 

properties (except for the binary diffusion 
coefficient) were taken to be constants, (II) 
the binary diffusion coefficient was taken to 
be inversely proportional to the local mixture 
density, (III) pressure gradient contribution 


to specie diffusion and momentum transport due 
to species diffusion were neglected, (IV) 
irreversible coupling effects between the tem- 
perature gradient and the specie concentration 
gradient (Soret and Dufour effects) were 
neglected, (V) gravity force was not considered, 
(VI) bulk viscosity was taken to be zero, (VII) 
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radiation heat transfer was not considered, and 
(VIII) chemical reactions did not occur. 
Furthermore, since the flow of interest is tur- 
bulent and present computing capabilities pro- 
hibit direct or large eddy simulations of 
turbulence, the conservation equations were 
density- weighted , ensemble averaged and closed 
by a K- € model of turbulence (28). The K-e 
model of turbulence employed was modified to 
account for some of the effects of compressibil- 
ity (29), streamline curvature (30), low 
Reynolds number (31), and preferential stress 
dissipation ( 32) . 

The dens ity- weighted , ensemble-averaged 
conservation equations and the K-e model of 
turbulence used in this numerical study are 
summarized in Table II. Equations (13) to (44) 
in Table II constitute a closed system in seven 
dependent variables: density (p), mass frac- 

tion of air (X^) , x-component of the velocity 
(u) , y- component of the velocity (v) , total 
energy ( e ) , turbulent kinetic energy (K) , and 
dissipation rate of turbulent kinetic energy (e). 

The boundary and initial conditions cor- 
responding to Eqs . (13) to (44) are given in 
Table III. The no-slip condition requires the 
fluid velocity next to a solid wall to be equal 
to the velocity of the solid wall. This condi- 
tion gives Eqs. (45) and (46). Here, it is 
noted that the K-e model of turbulence 
employed was modified to account for some of the 
effects of low-Reynolds number so that the model 
can be applied all the way to the wall (i.e., 
wall functions are not needed) . 

The gas temperature next to a solid wall 
equals the temperature of the solid wall. With 
all walls maintained at given temperatures, this 
boundary condition results in Eqs. (47) and (48). 

The walls are assumed to be nonporous. 

This boundary condition is expressed by 
Eq. (49). 

The turbulent kinetic energy and the dis- 
sipation rate of turbulent kinetic energy are 
both set equal to zero at solid walls because 
of the turbulence model employed. These two 
boundary conditions are given by Eq. (50). 

The stagnation temperature (T 0 ) and pres- 
sure (P Q ) of the air entering the combustion 
chamber through the intake port are constant as 
reflected by Eqs. (51) and (52). The mass 
fraction of the air entering the combustion 
chamber through the intake port is unity and is 
given by Eq. (53). At the intake port, the 
x-component of the velocity is zero and the 
y- component of the velocity is determined by the 
Bernoulli equation. These two boundary condi- 
tions are given by Eqs. (54) and (55). The 
turbulent kinetic energy and the dissipation 
rate of turbulent kinetic energy of the air 
entering through the intake port are given by 
Eqs . (56 ) and (57) . 

At the fuel injector, gaseous fuel (octane) 
is injected into the combustion chamber when the 
crank angle 0 is between 0^ and 0 2 . 

During fuel injection, the stagnation tempera- 


ture of the fuel (Tf) is kept constant as 
shown by Eq . (58). The fuel is injected at an 
angle a into the combustion chamber (see 
Fig. 1). The xand y-components of the velocity 
are given by Eqs. (59) and (60). The first- 
derivative of the pressure normal to the fuel 
injector is taken to be zero. This boundary 
condition is given by Eq. (61). The turbulent 
kinetic energy and the dissipation rate of tur- 
bulent kinetic energy of the fuel at the fuel 
injector are given by Eqs. (62) and (63). 

At the exhaust port, the back pressure 
(P e xh) i- s maintained constant. This condition 
is used to calculate the y-component of the 
velocity for the gas leaving the combustion 
chamber as indicated by Eq. (64). At the 
exhaust port, the x-component of the velocity 
along with the first-derivative with respect 
to y of the following variables are zero: 
mass fraction of air, pressure, temperature, 
turbulent kinetic energy, and dissipation rate 
of turbulent kinetic energy. These boundary 
conditions are given by Eqs. (65) and (66). 

At time t equal to zero, the rotor is at 
the beginning of the intake process (0 = 0) 
and the combustion chamber is filled with air 
at static pressure P^ and static temperature 
T^ . At this time, the air is in a state of 
homogeneous turbulence with zero mean velocity. 
The turbulent kinetic energy and the dissipation 
rate of turbulent kinetic energy of the air are 
both constant. These initial conditions are 
specified by Eqs. (67) to (72). 

NUMERICAL METHOD OF SOLUTION 

Solutions to the governing equations formu- 
lated in the above section can only be obtained 
by numercial methods. In this study, solutions 
were obtained by a finite different method. 

All finite-difference methods involve two 
approximations. First, the continuous domain 
of the problem must be replaced by time levels 
and a system of grid points. Second, the 
governing partial differential equations must 
be replaced by finite difference equations. 

These two approximations are described below. 

TIME LEVELS AND GRID SYSTEM - Finite- 
difference methods only provide solutions at the 
time levels and the grid points. The time 
increment between two successive time levels 
(i.e., the time-step size) must be made small 
enough to ensure numerical stability and tempo- 
ral accuracy. For the present problem, the 
rotor is started impulsively from 0 to 5000 rpm. 
For such a problem, it is necessary to use a 
very small time -step size in the beginning 
(Courant number much smaller than unity) . As 
time progresses, the time- step size used can be 
allowed to increase steadily until it reaches 
some preassigned value (corresponding to a 
Courant number much larger than unity) . Such 
variable time-step sizes were employed here. 

The time- step size used between time levels n 
and n+1 (denoted by At n ' 1 ' 1 ) was as follows: 
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TABLE II. - GOVERNING EQUATIONS 21 


Equation 


3 p 3 — ~ 3 — ~ _ 

£ + ^ pu + ^ pv = 0 

3 a 3 - — 3 a ~ 3 3 ~ 

^ * X A + ^ pUX A + ^ * VX A - « C — x ‘ + “ C - “ x 


3t 


3x 

3 


ay 


3X e 3x A 3y e 3y A 
3 t 

~ V T 4 ^ . 

3x e xx 3y e yx 

3_ t 3_ t 

ay '* ■ 3 T ax p e T xy + 3y v e T yy 


a a -~2 a — a - 2 - a t a t 

DU 4- pU 4- — pUV = - T (P 4 ~ pK) 4- T~ \i T 4- T u T 

“ 3y 3x 3 Ay yy A v 


3 3 - — 3 —2 a - 2 - _ 3 t 3 t 

at pv + ix pVU + ^ pv = ~ ^ (p + * pK) + ^ p ~ T — + ^ i*- T « 


ft + k (i + ? + 1 pK>:; + fc (i + ? + f pK) * 

3 , t ~ t - a , t ~ t~ 

3x v e xx xy 3y e yx yy 


+ jj-x 

3x e 


+ f-x f + |_ 

3y e 3y 3y 


pK + fc pSk + P " K 


f + f-|c |7h°- 
3x 3x | e l A 

*i -<)*(< 


+ ‘ C PA C PB 


T - T. 


*) 


ax *a 



C - c 
PA PBj 


*) ( ? - * 0 ] fe **} 

fc[R)*H[R)«] 


Eq. No. 

(13) 

(14) 

(15) 

(16) 


(17) 


3 - 


P€ 4- — puc + — PV€ = — y + 


3t ^ 3x 
where 


ay 


a_ \( Ht \ a* "1 a_ |Y \\ ii 
ax + o t j ax J + ay + « £ j ay 


~ 1 G _ P £ + D 

(18) 

— 1 I + s (19) 


= 2^-1? • V 
3x 3 ~ 

(20) 

= 2 f - f V • V 
ay 3 

(21) 

t 3u 3v 

yx _ 3y 3x 

(22) 

V = f + f* 
- 3x 3y 

(23) 

= y + P t 

(24) 

= C f pK 2 /c 
y y 

(25) 

= pD^ + c t 

(26) 

= V Sc t 

(27) 

= A + X fc 

(28) 
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TABLE II. 


Continued . 
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TABLE II. - Concluded. 


*6 = density-weighted mean of g, g - unweighted mean of g, p = den- 
sity, u = x-component velocity, v = y-component velocity, e = total energy, P 
« pressure, T = temperature, K = turbulent kinetic energy, c = dissipation rate 
of turbulent kinetic energy, y = dynamic viscosity, D^g = binary diffusion coef- 
ficient, X = thermal conductivity, Cpi = constant pressure specific heat of 
component i, = molecular weight of component i, hf = specific enthalpy 
of component i, * reference temperature , R - universal gas constant, Sc*. = 
turbulent Schmidt number = 1, Pr*. - turbulent Prandtl number = 1, = 0.09, 

C X = 1.44, C 2 = 1-92, <y K = 1.0, o c = 1.3, c[" = 2.24/C^ 3/3C = derivative 
normal to the boundary, fg^ equals 1 in the boundary layer next to solid walls 
and equals 0 outside the boundary layer, f^ = 1.0. 


TABLE III. - BOUNDARY AND INITIAL CONDITIONS a 


Location 

Equation 

Eq. No 

Housing surface 

u * V = 0 

(45) 


- dX 2 - 3Y 2 


Rotor surface 
and apices 

u = it - "’ v = 3t~ 

(46) 

Housing surface 

T = 300 K 

(47) 

Rotor surface 

T = 300 K 

(48) 

and apices 

3JL 


Housing and rotor 
surfaces, apices 

o 

it 

(49) 

Housing and rotor 
surf aces , ap ices 

K = € = o 

(50) 

Intake port 

P = l.oixio 5 N/Jti 2 
0 

(51) 

Intake port 

T = 300 K 
0 

(52) 

Intake port 

X A = 1.0 
A 

(53) 

Intake port 

u = o 

(54) 

Intake port 

? = - C DI [ 2(P o - ? >^] 1/2 

(55) 


/ „\2 


Intake port 

K = (yj 

(56) 

Intake port 

, - «•"/(.*) 

(57) 

Fuel injector 

T * 300 K 

(58) 

Fuel injector 

~ ( . e ' 6 i \ . 

(59) 

V o e 2 - e x ; 
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TABLE III. - Concluded. 


Fuel injector 

u = v/cotan a 





(60) 

Fuel injector 

= 0 
3C 





(61) 


_ 2 






Fuel injector 

K = (n 3 v) 





(62) 

Fuel injector 

c = K 3/2 ^n 4 L 2 j 

f / 

1 

\ 

1 1 

1/2 


(63) 

Exhaust port 

V - C DE [ 2 ( P - 

P exh) 

r J 



(64) 

Exhaust port 

U - 0 





(65) 

Exhaust port 

§- x « f- p = 

3y A 3y 

ay 

= f-K = 
3y 

f-* = 0 

3y 

(66) 


/; = p i/(RT./M A ) 


Everywhere inside 
the combustion 
chamber at t = 0 


X A = 1 


u = v * 0 




e = P 

K = constant 
c = constant 


(67) 

( 68 ) 
(69) 


i[( M A C PA- R )j R+ ( h A- C PA T i)] (70) 


(71) 

(72) 


= 0.03, n 2 - 0.0075, = intake port diameter, V 0 = 40 m/sec, n 3 

^ 0.03, n A = 0.0075, L 2 = fuel injector diameter, C DI = discharge coefficient 
for intake port, C D p = discharge coefficient for the exhaust port, Pf = 10 5 
N/m 2 , P exh = 9.9x10 s N/m 2 . 


[9(n + 1) - 6(n)]/Q n = 0,1,2 n 


ntl 


At 


where 


s 

(73) 


[A0 f /ft n = n g + 1, n g + 2, . . . ,n t 


2 3 

6(n) = 0 + a x (n/n g ) + a 2 (n/n g ) + a 3 (n/n g ) 


a* = A0 + A0_ - 2 (0 r - 0 )/n / 1/n 

3 f 1 t o sj 




(74) 


+3/n 2 + 2/n 3 (75) 

s s 


[“l - 


(6 - 6 )/n - a (1/n 

r o s 3 s 


1/n ' ) ] / [ 1/n ® ” 1/n s (76) 


A0. 


3 1 = 9 f - ®o - a 2 - a 3 

■ [‘, - <e f * V] {“t - ",] 


(77) 

(78) 


In Eqs . (73) to (78), ft is the angular speed 
of the shaft (ft = 5000 rpm) . 0 O is the crank 
angle at which we begin computations (0 O - 0) . 
0f is the crank angle after which the time- 
step size becomes a constant (Of = 1 radian). 
The parameter n s denotes the time level when 
0 = Of (n s = 800). A8i is the crank angle 
traversed by the shaft during the first time 
step (A6x = 0.001 A0f). A0f is the crank 
angle traversed by the shaft during each time 
step after time level n = n s or crank angle 
0f. The parameter n^ is the total 
number of time steps needed to traverse 6ir 
crank angles (n^ = 10 000) . 

For computational efficiency, the number 
of grid points (intersections of grid lines) 
used should be kept to the minimum required to 
resolve spatially all significant features of 
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the flow. Here the system of grid points 
employed was as follows. In the x-y-t or 
Cartesian coordinate system (Fig. 2(a)), the 
grid points moved with the combustion chamber 
as the rotor rotated inside the housing. This 
reduced the number of grid points needed. To 
further reduce the number of grid points needed, 
grid points were clustered in regions where they 
are needed and widely spaced in regions where 
they are not needed. Here, more grid points 
were placed near the rotor, housing, and seals 
to resolve the steep gradients of the dependent 
variables in those regions. Finally, grid lines 
interesected the rotor and the housing perpen- 
dicularly to facilitate implementation of normal 
derivative boundary conditions. 

To facilitate derivation of finite- 
difference equations, implementation of boundary 
conditions, and programming, the nonuniformly 
distributed and moving grid points in the x-y-t 
coordinate system were mapped onto another 
coordinate system, the £-n-T coordinate system. 
In addition, grid lines in the x-y-t coordinate 
system corresponded to coordinate lines (i.e., 
lines of constant £ or constant n) in the 
coordinate system. 

This mapping of grid points is known as 
grid generation. Here, the mapping or grid 
generation was achieved by transfinite inter- 
polation (33,34). The functional relationship 
between the two coordinate systems obtained by 
transfinite interpolation was derived in 
Ref. 34 and summarized by Eqs. (79) to (89) in 
Table IV. 

FINITE-DIFFERENCE EQUATIONS - The finite 
difference equations used to approximate the 
governing equations summarized in Table II were 
derived by the implicit-f actored method of Beam 
and Warming (35,36). The details of the finite- 
difference equations used along with their der- 
ivation were described elsewhere (37). Here, 
only the essence will be presented. 

The first step in deriving the finite- 
difference equations was to replace all deriva- 
tives with respect to t, x, and y by derivatives 
with respect to t, £, and r\. The second 
step was to replace all derivatives with respect 
to x (including the metric coefficients from 
the coordinate transformation) by the second- 
order accurate in time Crank- Nico Ison time- 
difference formula. The third step was to 
employ a linearization procedure. The fourth 
step was to employ an approximate factorization 
procedure followed by splitting of the differen- 
tial operator. Finally, all spatial derivatives 
(including the metric coefficients) were 
replaced by second order accurate in space 
central-difference formulas. Here it is noted 
that all metric coefficients were evaluated by 
finite difference formulas instead of by 
analytical differentiation. This was done to 
avoid geometrically induced errors (34,38). 

The finite- difference equations derived by 
the steps described above formed systems of 
linear equations with block tridiagonal 
coefficient matrices and were solved by the 
noniterative Thomas' algorithm. These finite- 


difference equations provided the solutions at 
every grid point and time level presented in the 
next section. 

RESULTS 

Numerical solutions were generated to 
investigate the flow patterns inside the com- 
bustion chambers of the rotary engine depicted 
in Fig. 1 during intake, compression, gaseous 
fuel injection, expansion, and exhaust. Den- 
sity, mass fractions of fuel and air, x- and 
y-components of the velocity, energy, tempera- 
ture, pressure, turbulent kinetic energy, and 
dissipation rate of turbulent kinetic energy 
were calculated at 1071 (i.e., 51-by-21) grid 
points (see Fig. 3) for 10 000 time levels. 

Only results for the velocities (velocities 
relative to the stationary housing) are pre- 
sented. These results are shown in Figs. 4 to 
16 illustrating the flow patterns. Do not infer 
velocity gradients from these figures because 
the velocity vectors were plotted to best 
illustrate the flow patterns. 

Figures 4 to 6 illustrate the flow patterns 
during the intake process. During the early 
part of the intake process, both the intake and 
the exhaust ports were open. Accordingly, air 
was found to enter the combustion chamber 
through the intake port and to exit through the 
exhaust port as shown in Fig. 4. Here, the 
stagnation pressure of the air entering the 
combustion chamber P Q was made slightly 
higher than the exhaust port back pressure 
p exh slightly turbocharged) so that 

back flow through the intake port did not qccur. 
Later in the intake process , only the intake 
port was open and a high speed jet issued into 
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TABLE IV. - COORDINATE TRANSFORMATION® 


Equation Eq. No. 

3* 2 (5,t) 

x(5,ti,t) = X^e.rjh^n') + x 1 (Z,-')h 2 (n’> - K 2 (5) — h^n 1 ) 

9Y (5,t) 

- K x (5) V*i’> (79) 

3X (5,t) 

y( 5 i ti , t ) = y 2 (5,T)h 1 (n') + 5^(5, -r)h 2 +■ k 2 (5) — ^ V v > 

3X (5,t) 

+ K(5)— h.(n') (80) 

1 35 4 


t = T 


(81) 


where 


h 1 (n’ ) 


2(n ' ) 3 - 3(n’) 2 + 1 


(82) 


h 2 (n’) = -2(n*) 3 + 3(n’) 2 


h 3 (n-) = (n ' ) 3 - 2 (n ’) 2 + n’ 


h.(n’) = (n') 3 - (n’) 2 

4 


(83) 

(84) 

(85) 


>• • 0 5 ( B , * ') - *, | l * (*, * jl(\ ~ ')] ""''’l 


(v 0 ‘ M‘ * [(v l )/(*. - 


(25-1) 


( 86 ) 
-1 


i -wz-° s ( b ^ i )* e 4 i+ [S* i )I( b ^ 1 )] <25 ' i> 1 


(87) 
-1 


1^(5) = k 2 (5) = 2 


|^x 1 (5.t) - x 2 (5,x)j 2 +|^y 1 (5,t) - Y 2 (^,T)j 2 | 


( 88 ) 
1/2 


(89) 


ax l» Y l» x 2> and Y 2 are Siven by Eqs. (1) to (12) in Table I with 
the parameters A and V related to 5 as shown by Eqs. (87) and (88) in this 
table; = 0.1; = 0.3. 


the combustion chamber (Figs. 5 and 6). In 
Fig. 6, a very small region of recirculating 
flow can be observed near the intake port. 
However, this recirculating flow region was 
found to appear aperiodically and then disap- 
pear* The disappearance of the recirculating 
flow is due to the general motion of the entire 
combustion chamber in one direction. 

Figures 7 and 8 i llustrate the flow pat- 
terns during compression. Recirculating flows 


were not observed during the compression 
process . 

Figures 9 to 11 illustrate the flow pat- 
terns during gaseous fuel injection. Here 
gaseous fuel (octane) was injected into the com- 
bustion chamber between crank angles 0^ = 7.4 
rad or 425° and 02 * 8.5 rad or 487°. The 
velocity of the fuel jet was given by Eqs. (59) 
and (60) in Table III. The maximum speed of the 
fuel jet was V 0 = 60 m/sec and the fuel was 
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Figure 9. - Flow pattern at crank angle 
0 - 7.59 radians or 435°. 




Figure 11. - Flow pattern at crank angle 
9 = 8. 27 radians or 474°. 


injected at an angle of a = 45° (see Fig. 1). 
When the fuel was first injected, its momentum 
was low in comparison to the momentum of the gas 
already inside the combustion chamber. Conse- 
quently, at this time the fuel jet could not 
penetrate very deeply into the combustion cham- 
ber as shown in Fig. 9. Later as the momentum 
of the fuel jet increased, it penetrated much 
deeper into the combustion chamber and a rel- 
atively large recirculating flow region can be 
observed as shown in Fig. 10. Towards the end 
of the fuel injection process, the momentum of 
the fuel jet was again low in comparison to the 
momentum of the gas inside the combustion cham- 
ber and the fuel jet was deflected as shown in 
Fig. 11. 
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Results for the fueL injection process 
should be interpreted carefully. First, in 
realistic rotary engines, a liquid fuel instead 
of a gaseous fuel is injected. A liquid fuel 
jet has much more momentum than a gaseous fuel 
jet at identical jet velocities. Second and 
more important, in realistic rotary engines, the 
fuel jet invokes a three-dimensional phenomenon 
known as a jet in a crossflow. Jets in cross- 
flows cannot be accounted for by two-dimensional 
numerical studies. 

Figures 12 and 13 illustrate the flow pat- 
terns during the expansion process. No recircu- 
lating flow patterns were observed. 

Figures 14 to 16 illustrate the flow pat- 
terns during exhaust. Here, we can observe some 
recirculating flows near the exhaust port as the 
gas rushed out of the combustion chamber. 

The results presented above were not com- 
pared with other numerical or experimental data 
because they do not yet exist. As noted in the 
Introduction, the present study constitutes the 
first attempt at numercially studying the 
unsteady, multidimensional flow inside rotary 
engine combustion chambers. 



Figure 15. - Flow pattern at crank angle 
0= 17. 1 radians or 980°. 




Figure 12. - Flow pattern at crank angle 
0= 11. 10 radians or 436°. 




Figure 14. - Flow pattern at crank angle 
0 = 16. 72 radians or 958°. 


Figure 16. - Flow pattern at crank angle 
0= 18.1 radians or 1037°. 

CONCLUDING REMARKS 

The results generated in this study illus- 
trate some of the major features of the flow 
patterns taking place inside the combustion 
chambers of a two-dimensional rotary engine with 
gaseous fuel injection. From the results 
obtained so far, it appears that the mathemati- 
cal model developed here can be used to study 
the effects of a number of engine design and 
operating parameters (such as those relating to 
the fuel injector, rotor pocket, engine speed, 
and in take -exhaust manifold pressure ratio) on 
the flow patterns. 

Clearly, a two-dimensional study cannot 
account for all the important factors which 
influence combustion inside a real rotary 
engine. In particular, the effects of a jet in 
a crossflow (taking place during intake and fuel 
injection) as well as Taylor vortices can only 
be accounted for by a three-dimensional analy- 
sis. Since jets in crossflows and Taylor vor- 
tices can significantly affect the flowfield 
inside rotary engines, a three-dimensional 
numerical study is already underway. 
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